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"o ! Abstract 



We propose a new, discretized model for the study of 3+1-dimensional canonical quan- 
tum gravity, based on the classical SL(2, (D)-connection formulation. The discretization takes 
place on a topological iV 3 -lattice with periodic boundary conditions. All operators and wave 
functions are constructed from one-dimensional link variables, which are regarded as the fun- 
damental building blocks of the theory. The kinematical Hilbert space is spanned by poly- 
nomials of certain Wilson loops on the lattice and is manifestly gauge- and diffeomorphism- 
invariant. The discretized quantum Hamiltonian H maps this space into itself. We find a 
large sector of solutions to the discretized Wheeler-DeWitt equation Hip = 0, which are la- 
belled by single and multiple Polyakov loops. These states have a finite norm with respect to 
a natural scalar product on the space of holomorphic SL(2, (D)-Wilson loops. We also inves- 
tigate the existence of further solutions for the case of the l 3 -lattice. - Our results provide 
for the first time a rigorous, regularized framework for studying non-perturbative quantum 
gravity. 
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1 Introduction 



All attempts to define a lattice discretization of four-dimensional quantum gravity, in 
analogy with the rather powerful methods employed in quantum chromodynamics, have 
been riddled with difficulties. A major problem is that of incorporating the diffeomorphism- 
invariance of the theory. The discretization typically destroys this symmetry already at the 
classical level, similar to the way in which a lattice discretization of Yang-Mills theory breaks 
its rotational symmetry. However, in the case of gravity this is much more serious, since the 
diffeomorphism group is infinite-dimensional, and acts non-linearly on the underlying mani- 
fold. A number of questions arise: (i) does the discretized quantum theory possess a residual 
diffeomorphism symmetry? (ii) how is the continuum limit to be taken? and (hi) what 
does the resulting theory look like? (This should be a non-perturbative, diffeomorphism- 
invariant description for quantum gravity in an "unbroken" phase.) The answer to (i) is 
usually "no" , whereas the other two questions are hard to address and have not yet found a 
satisfactory solution within either Lagrangian or Hamiltonian approaches. Using Euclidean 
path integral methods, a main problem is that of finding the correct measure, without over- 
or undercounting the physical configurations (currently the most active research programs 
are quantum Regge calculus and dynamically triangulated gravity; see [1] for a recent re- 
view). A central problem in the canonical approaches is that of finding a closing algebra of 
diffeomorphism constraints on the phase space of the regularized theory (see, for example, 
the discussion in [2]). The unresolved status of the closure problem is the reason for which 
many Hamiltonian discretized models have not advanced very far [3,2]. 

The starting point of this paper is a particular form of Hamiltonian lattice gravity. How- 
ever, instead of implementing the action of the (three-dimensional) spatial diffeomorphisms 
directly on the phase space of the theory, we will use a framework in which diffeomorphism 
invariance is manifest. This approach is inspired by the loop space formulation of canonical 
quantum gravity, based on the classical reformulation of Einstein gravity in terms of a pair 
{A l a ,Ef) of SL(2, (D)-Yang-Mills variables due to Ashtekar [4,5]. In the original paper by 
Rovelli and Smolin [6], the wave functions ^(7) of the quantum theory are labelled by spatial 
loops 7, and the diffeomorphism invariance is formally implemented by selecting those wave 
functions that are invariant when the loop argument 7 is moved by a diffeomorphism. Physi- 
cal wave functions therefore depend on diffeomorphism equivalence classes [7] of closed loops. 
Using these ideas, solutions to all of the quantum constraints, including the Wheeler-DeWitt 
equation Hip = 0, have been found [6] (for an overview of this and other solutions, see [7]). 
This involves the choice of a regularization and factor-ordering for the quantum Hamiltonian 
constraint H, and is rather formal, in the sense that there is no well-defined scalar product 
on the space of solutions, and little control over the influence of different regularizations and 
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factor orderings on the structure of this space. As a result the status of these solutions has 
remained unclear and the issue of reality conditions could not be addressed. In the Ashtekar 
formulation in terms of complex canonical variable pairs [A l a , Ef ) , such conditions have to 
be implemented in the quantum theory to make sure that real, and not complex gravity is 
described. 

Recently, there have been proposals for defining continuum loop representations rigor- 
ously , regarding them as non-linear analogues of the Fock representation based on quantum 
loop states (see [8] for a review and further references). Central to this line of research are 
the construction of suitable domains for wave functions depending on connections modulo 
gauge transformations and diffeomorphism-invariant measures on such spaces. However, so 
far these efforts have addressed only the kinematical structure of the quantum theory, without 
incorporating any of the dynamical issues regarding the quantum Hamiltonian. A rigorous 
analysis of the structure of the solution space to the Wheeler-DeWitt equation and of phys- 
ical operators on this space is therefore still lacking. In the present work we will suggest a 
discretized version of canonical quantum gravity that can deal with these issues. 

The discretization will take place on a periodic cubic lattice, but this is not to be thought 
of as a fixed lattice embedded in physical, Euclidean three-space, but merely as a topolog- 
ical quantity, defined by those of its properties that would remain unaffected by a smooth 
diffeomorphism. It may thus be thought of as representing a diffeomorphism equivalence 
class of cubic lattices. All physical quantities are defined in terms of link variables, and 
the discretized Hamiltonian acts purely combinatorially on wave functions labelled by lattice 
loops. The formalism is manifestly diffeomorphism- and therefore also scale-invariant. The 
basic link variables of our formulation are a Kogut-Susskind pair of an SL(2, (D)-valued link 
holonomy and a corresponding sl(2, (D)-momentum variable. However, the corresponding ba- 
sic operators in the quantum representation are not self-adjoint. This is acceptable, since 
they do not correspond to any physical observables. It comes about since we propose to 
implement the reality conditions, following the ideas of [5,3], as a holomorphicity condition 
on wave functions (this leads to the correct number of degrees of freedom in the quantum 
theory). Also the Hamiltonian operator turns out to be non-selfadjoint. This is no real reason 
for concern either, if one accepts the argument that the self-adjointness of the Hamiltonian 
in a generally covariant theory is not a strict physical requirement (for related discussions, 
see [9]). 

On a cubic iV 3 -lattice with periodic boundary conditions, given a specific discretization 
of the Hamiltonian H together with a particular operator ordering for H, we are able to find 
an infinite set of solutions to the Wheeler-DeWitt equation, that moreover have finite norm 
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with respect to the inner product induced from the original kinematical Hilbert space. They 
are labelled by so-called Polyakov loops and multiple Polyakov loops, which are well-known 
from their role as order parameters for the phase structure of lattice gauge theory (see, 
for example, [10]). Their appearance in the context of lattice gravity is curious, although 
they play here a quite different role, namely, that of parametrizing the solution space to the 
Wheeler-DeWitt equation. These solutions are the lattice analogues of the non-intersecting, 
smooth- loop solutions found by Rovelli and Smolin [6]. This for the first time provides a 
regularized model for non-perturbative quantum gravity that can be used to study physical 
observables. 

The paper is organized as follows. In the next section we introduce all the necessary 
ingredients for setting up the Hamiltonian lattice theory, and explain some features of the 
holomorphic representation we will be using for constructing a scalar product on SL(2, (D)- 
wave functions. We derive an important explicit formula for relating arbitrary polynomials 
in the SU(2)- and SL(2, (D)-theories. In section 3, we establish a relation between the zero- 
eigenvalue solutions of our quantum Hamiltonian and another Hamiltonian, induced from the 
SU (2)-theory. This suggests a close connection between the real and the complex, holomor- 
phic theories. Next we demonstrate that there is a large subspace of the Hilbert space that is 
annihilated by the Hamiltonian constraint. In section 4, we take a closer look at the special 
case of the l 3 -lattice and explore the possibility of finding solutions beyond the ones labelled 
by the Polyakov loops. We illustrate some of the technicalities involved, without being able 
to identify any solutions explicitly. The last section contains our conclusions and an outlook. 



2 The general formalism 

We first recall the basic Hamiltonian variables for the S'[/(2)-lattice gauge theory, before 
discussing the complexified framework for SL(2, C) = SU{2)<^. This is necessary for setting 
up a discretized version of the connection formulation of canonical gravity. - With each lattice 
link we associate an element g G SU(2), parametrized by a matrix V(g) in the defining two- 
dimensional representation as 

y b _ ( cto + icti ct2 + ia3 
\ — «2 + icx$ ceo — ia\ 

with oti £ IR, and subject to the condition X]i=o a ? = 1- ^he ma t rrx Va B can be written as 
a (real) linear combination of the unit matrix 11 and the three r-matrices defined by 
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n = (j °). » = (_», J), -, = (» J). (-) 



The r-matrices satisfy [t^t,-] = 2eijkTk- The differential operator 



/ 5 n s _ 1 ( d - id 1 -d 2 -id 3 \ 



acts on the representation matrices (2.1) and satisfies 



{-^)a B Vc D = \8 a d 8c B . (2.4) 

The operators corresponding to the classical momentum variables pi (with a gauge algebra 
index i) are given by 



k = - iT%A B VB C { JL )c A (25) 



and satisfy 



(2.6) 



h{v- i ) A c =;(nA, B 



In terms of coordinates, we have 



2, 

P2 = -^(a 2 d - a 3 di - a d 2 + aid 3 ) (2.7) 
P3 = ^(^do + ol>$\ - a x d 2 - a d 3 ). 

Defining the operator V as multiplication by the matrix V therefore leads to the commutation 
relations 
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[Va B ,Vc D ] = 

[P^V A c ] = - l -n A B v B c (2.8) 
\pi,Pj] = ie ijk p k , 

which are the quantum equivalents of the Poisson brackets of the corresponding classical 
quantities). These commutation relations are familiar from the Hamiltonian SU (2)-\&ttice 
gauge theory [11], where the pi are hermitian operators. Here, in contrast, we shall associate 
three complex (i.e. SL(2, C)) degrees of freedom with each link of the hypercubic lattice. As 
usual, operators (V,p) associated with different links commute. 

The entire construction (2.1-8) makes sense also if we complexify the group to SU (2)@ = 
SX(2,C). We will denote the corresponding group parameters by complex numbers af, 
i = . . . 3, again subject to a condition Yl^=o( a f) 2 = 1- The operators V and p are taken to 
act on a space of holomorphic wave functions on the group manifold of SL(2, <E), whose inner 
product will be specified below. Note that due to the non-compactness and non-abelianness 
of the group SL{2, (D), there is no analogue of the bi-invariant Haar measure dg, which exists 
on SU (2). However, thanks to the work of Hall [12], we know there exist unitary holomorphic 
transforms from the space L 2 (SU(2),dg) of square-integrable functions on SU(2) to spaces 
L 2 (SL(2, (D), dv)^ of holomorphic and i/-square-integrable functions on SL(2, C). (These 
are analogous to the Segal-Bargmann integral transform from L 2 (IR n ) into the holomorphic 
functions on <C n .) This provides us with the desired scalar product on functions on a complex 
group manifold. Note that Hall's results have recently been used to construct a coherent state 
transform for spaces of connections in the continuum [13]. 

We now recall some details of Hall's construction insofar as they are relevant to the present 
discussion. For each real t > 0, there is a "coherent-state transform" C t ■ L 2 (SU (2) , dg) — > 
L 2 (SL(2,(E),d» t ) n defined by 

[C t (f)](gv) ■= [ f(g)Pt(g- 1 g< c )dg, (2.9) 

J SU{2) 

where / G L 2 (SU(2), dg), g G SU(2), g^ £ SL(2, C), and p t is the heat kernel for the Casimir 
operator A = —^YliPi on SU(2), i.e. the fundamental solution at the identity of the heat 
equation dp/dt = ^Ap. More precisely, since the argument of pt in (2.9) is a complex group 
element, we are using its analytic continuation, which is well-defined (see [12] for details). In 
terms of the explicit parametrization (2.1) for the matrices V(g) and the normalized Haar 
measure dg, and using a series expansion for the heat kernel, one obtains 
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i II 

(2J + 1) e-^' +1 )*/ 2 U 2j (a a$ + ai af + a 2 a% + a 3 a^) = 

l \f l ~ a l v /l_a o- a ? 
da / dai / afa 2 V (2j + 1) 

1 -y/\-a%-a\ '-°>2- 
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(2.10) 



(whenever the infinite sum over j on the right-hand side converges), where the U 2 j denote 
the Chebyshev polynomials of the second kind. The image of a square-integrable function 
f(g) on 577(2) is a holomorphic function on 5L(2,(D), square-integrable with respect to du t , 
which is essentially the heat kernel measure on the quotient space SL(2, C)/5C7(2). 

Since our aim is a manifestly gauge-invariant description of lattice gravity in a holomor- 
phic loop representation, several issues have to be resolved. Firstly, we are not aware of a 
simple explicit expression for the measure dv t in the holomorphic representation obtained via 
(2.9), and therefore have to look for functions f(g) with a simple transformation behaviour, 
preferably elements of an orthonormal basis of L 2 (SU (2) , dg) , which will be mapped into or- 
thonormal functions of L 2 (SL(2, (D), dv t )^, since C t preserves scalar products. Secondly, we 
will work with functions that are gauge scalars and can be expressed as functions of Wilson 
loops, i.e. traced holonomies of closed loops on the lattice. The lattice Hamiltonian maps 
such functions into themselves. 



To obtain the holomorphic transform for general 5C/(2)-wave functions on the lattice, we 
need to take the product over all lattice links of the transform for a single link, formula (2.10). 
First however we will give an (overcomplete) set of functions of a single copy of 577(2), i.e. 
on a single link, that have a simple transformation behaviour under the transform (2.10). 
They are given by appropriate sums of polynomials in the four real parameters CKj (restricted 
to the submanifold 577(2) C M 4 ) and can be labelled by the exponents in the polynomial 
ctQ° a™ 1 (*2 2 a^ 3 of highest order occurring in the sum. One finds 
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p(n ,n 1 ,n 2 ,n 3 ) := 




- 2 2 Oo-i)( JO -l)!(n -2(jo-l))! 




(n )!(n - Jo + ni + n 2 + n 3 + 1)! 



(n + n\ + n 2 + n 3 )! 




, 2 20 1 -i)(j 1 -l)!(n 1 -2( Jl -l))! 




(rai)!(n - jo + ni - ji + n 2 + n 3 + 2)! 



(«o - Jo + ni + n 2 + n 3 + 1)! 




(n2)!(n - jo + "-1 - ji + n 2 - j 2 + n 3 + 3)! 



(n 3 )!(n - jo + ni - ji + n 2 - j 2 + n 3 - j 3 + 4)! 



(«o -jo + ni -ji + n 2 -j 2 + n 3 + 3)! 

(2.11) 



(no ~ jo + ni - Ji + n 2 + n 3 + 2)! 



For the holomorphic transform of p(no, n±, n 2 , n 3 ) one obtains 



[C t (p(no,n 1 ,n 2 ,n 3 ))](af) = e -(«O+ni+n 2 +n3)(n +n 1+ n 2 +n3+2) t /8 p(n0) ^ „ 2> ng) C ) (2J2) 



where by p(no, ni, n 2 , n 3 ) we denote the expression (2.11) with the real parameters a, re- 
placed by the corresponding complex quantities af . That is, to find the (inverse) holomorphic 
transform of a polynomial function of the on (af), one first has to express it as a linear combi- 
nation of the p(no, n\, n 2 , n 3 ) (p(n , n\, n 2 , n 3 ) c ) and then use (2.12). Both (2.11) and (2.12) 
are crucial formulas for relating the SU(2)- and the holomorphic SL(2, C)-representation, and 
will be used in the following sections. The next step is the identification of gauge-invariant 
combinations of the p(no, ni, n 2 , n 3 ). If the configuration space consisted of just one link 
(with endpoints identified), our task would be straightforward: all gauge-invariant quantities 
one can construct in that case are functions of ^Tr V = ao, and a complete basis is given by 
the Chebyshev polynomials ?7 no (o;o)- 

The situation on the hypercubic lattice is more complicated since general gauge-invariant 
quantities are functions of traces of holonomies around arbitrary lattice loops, which do not 
necessarily factorize into products of link contributions. Moreover, there is the additional 
problem of finding a set of gauge-invariant functions on XiSU(2) that is complete (i.e. spans 
the Hilbert space of square- integrable functions), but at the same time not overcomplete, i.e. 
free of redundant degrees of freedom. This is a well-known complication with intrinsically 
gauge-invariant formulations of lattice gauge theory, and there are various strategies of dealing 
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with it. In the section 4 we will address some of these difficulties in the context of the lxlxl- 
lattice. 

Formula (2.11) can be used to construct gauge-invariant functions with a simple transfor- 
mation behaviour that are labelled by lattice loops. This is important because the (overcom- 
plete) basis of gauge-invariant functions {TV V^V^ ...Vi n , 7 = l^ol^o. . .ol in a lattice loop} 
often appears in applications. Using the parametrization (2.1) for the link matrices Vi, 
Tr Vi-J/iz . . . Vi n is a homogeneous sum of polynomials in those parameters. Re-expressing in 
each summand the link contributions in terms of the functions p(uq, n\, 112, 713), one obtains 
quantities that transforms like (2.12) with t-dependent exponential factors. 



3 Solutions to the Wheeler-DeWitt equation 

Next we study the action of a discretized form of the phase space Hamiltonian for canon- 
ical gravity. Let us label lattice sites by an integer n and the three positive directions 
emanating from each site by a = 1, 2 1 or 3. Thus the canonical variables are given by V(n, a) 
and p(n, a). We denote the holonomy of a plaquette loop based at the site n in the d-6-plane 
by V(n, P. g ), that is, V(n, P ah ) = V(n, a)V{n + a, b)V(n + b, a)- l V(n, S) -1 . 

We require that in the limit as all link lengths go to zero the continuum Hamiltonian is 
recovered. Since our lattice was assumed to be purely topological, we define this limit with 
respect to an auxiliary Euclidean coordinate system (the three lattice directions coinciding 
with the three coordinate axes) in which all links have length a. As a — > 0, one derives the 
usual expansion for the plaquette holonomy 

V(n, P hh )A B a -^° t A B + a 2 F k ah r kA B + 0(a 3 ), (3.1) 

where F a b is the a-6-component of the field strength associated with the selfdual connection 
A. For the momentum variable Pi(n, a) we require that 

Pi {n, a) ^ a 2 E?(n) + 0(a 3 ) (3.2) 

in terms of the continuum momentum density E, since p is like a momentum variable smeared 
in one (out of three) spatial directions. Therefore, if we associate with each lattice site n the 
lattice Hamiltonian 
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# c (n) = Y,e ljk Mn,a) Pj (n,b)TT(V(n,P. i )r k ), (3.3) 

d<b 

this to lowest order leads, up to a power of a, to the correct continuum limit 

H c ( n ) a & ^ k ~EfE b 3 F abk + 0(a 7 ). (3.4) 

The total Hamiltonian is given by the sum ^2 n H^{n). As a result of the "dimensional 
transmutation", the discretized Hamiltonian is a topological quantity, in the sense that it 
depends only on link variables, with no reference to the link length a, and acts on loop wave 
functions in a combinatorial way. As long as we do not re-introduce a length scale, the 
formulation is therefore purely topological. 

Just as in lattice gauge theory, the requirement of the correct continuum limit does not fix 
the discretized Hamiltonian uniquely. In the present work we do not investigate the question 
whether a different choice of H®(n) leads to equivalent results. In going to the quantum 
theory, another ambiguity arises in the choice of the operator ordering of ij c (n). The most 
commonly used operator ordering for the Hamiltonian is the one with the operators pi to the 
right, but also the opposite ordering with both of the pi to the left is sometimes used [14]. 
In our present investigation we will be using the former, i.e. 

k<C = E E eljfcTr ^(n, P &i )r k ) Pi(n, a) p j (n, b). (3.5) 

n a<b 

Since the spatial diffeomorphisms have been taken care of, the only remaining task is to 
look for holomorphic wave functions ^ c G X;L 2 (SX(2, C), dv t )^- that solve the analogue of 
the Wheeler-DeWitt equation, 

H € i) € = 0. (3.6) 

The existence of solutions depends on the spectral properties of the Hamiltonian oper- 
ator . Since the constituent operators V and p are not selfadjoint in the holomorphic 
representation, we do not have any a-priori information about the spectrum of H . There 
are three possibilities: 
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(i) = has non-trivial solutions in X/L 2 (SX(2, ©), dv t ) n ; 

(ii) H^tp® = does not have solutions in X; L 2 (SL{2,V),dv t ) n , but can be solved "as a 
differential equation", i.e. there are non-square-integrable solutions; 

(hi) there are no solutions. 

In case (i), the solution space inherits a scalar product from the original Hilbert space, whereas 
in case (ii) one still has to define a suitable inner product. - A number of identities are useful 
in computing the action of the Hamiltonian on gauge-invariant wave functions (which all 
contain terms of the form Tr V^V^ . . .). The first one is 



Jjk^ B D D r B B s: D in 7 \ 

£ TjA T k C =TiA OC —Tie $A , (3.7) 

from which follow two identities for the traces of holonomies, 

e^Tr (VcTjVpn) = T> VpTr (V a n) - Tr V a Tx (Vpn) 
e^' fc Tr (V a Tj)Tr {V P T k ) = Tr {V a nV p ) - Tr (V^n), 

where a, (5 are two lattice loops intersecting at n, and the product loop a o (3 is obtained by 
the usual loop composition at n. Secondly, we have 



r lA B T lC D = 5 A B 5c D -2S a D S b c . (3.9) 

Lastly, there is the well-known identity relating a product of two Wilson loops of SX(2,C)- 
holonomies to a sum of two Wilson loops, 



Tr y Q Tr Vp = Tr V a0 p + Tr V ao/3 -i = Tr V a V p + Tr V a V^\ (3.10) 

By virtue of these identities, the action of the Hamiltonian on states of the form Tr V^V^ . . . 
may be interpreted as cutting and joining of the lattice loop arguments, as is typical for 
the loop representation. To what extent such a geometric interpretation is useful in finding 
solutions to the zero-eigenvalue equation depends to some extent on the type of basis chosen 
for the quantum states. As in the case of Hamiltonian lattice gauge theory, the Hamilto- 
nian couples neighbouring links due to the occurrence of the plaquette holonomy operators 
V>,R g ) in H{n). 
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Next we show that solutions to the zero-eigenvalue equation H^ip® = are in one-to-one 
correspondence with solutions to Hmdip^ 1 = 0, where Hi n d is the (self-adjoint) Hamiltonian 
induced from the £f/(2)-lattice theory. That is, we take as in (3.5), and substitute the 
operators by their real counterparts acting on XiL 2 (SU(2), dg) (i.e. take all parameters ctj 
etc. to be real), before translating it to the holomorphic representation using the transform 
(2.10). 

Let us for the moment assume we are given a complete orthogonal basis of wave functions 
{x(n), n = (ni,n2, ■ ■ ■ nd)}, d = dim(x/SI/(2)/ x s SU(2)), for the gauge- invariant (real) 
Hilbert space L 2 (xiSU(2)/ x s SU(2), 7r(]X dg)), where tt(Y\i dg) denotes the projection to 
the quotient space of the product of Haar measures, I is the number of links and s the number 
of sites. Its elements are labelled by integers nj and transform according to 



Ct(x(n)) = e-^ X C (n), (3.11) 

and ft(n) > 0. Since the exponential factors are just real numbers rescaling the basis 
(which we did not assume to be orthonormal), we can write any square- integrable ele- 
ment of the holomorphic Hilbert space as a real linear combination of the X C (^)> an d have 
< (fn) j X^ ip) >~ Sm,n- The action of the Hamiltonian H acting on a general vector 
J2n a (^0 X (^)j a ("0 ^ IR> can be written as a matrix equation 



n ft fn 

The Wheeler-DeWitt equation (3.6) is therefore equivalent to an infinite tower of equations 
for the coefficients a(n), 



^M(n,m)a(n) = 0, Vm. (3.13) 



The norm of a vector Yln a (™)X (™) can be calculated using its inverse image under 
(3.11), and is given by 
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By contrast, the Hamiltonian Hf nd acts on holomorphic wave functions according to 

H? nd E b(n) X C (n) = E E e " /t( " ) M ("' ™) eMH) b ^ * V)> ( 3 - 15 ) 

n rim 

leading to the set of conditions 

J2M{n,m)e ft ^h(n) = 0, Vm (3.16) 

n 

on the coefficients b(ft) of zero-eigenvectors. Clearly they are related to the solutions of (3.13) 
by a(n) = e^ t(jl ^b{n). However, the norm of the solution vector is in general different, and we 
have 



iiE 6 ^)x c ^)ii^(2,c) = JEw^^iix^iiLp) < iiE a (")x c («)iUL(2,o- 

n V n n 

(3.17) 

It may therefore happen that a zero-eigenvalue solution to J2n a (^) X C (^) = is square- 
integrable, while the corresponding solution to Yin a(n) x C ("0 = is not. 

The form of the quantum Hamiltonian (3.5), like its counterpart in the continuum the- 
ory, is sufficiently complicated so as not to lead us to expect the eigenvalue problem could 
be soluble trivially. What comes to help in the analogous problem in the continuum is 
the existence of "algebraically special" solutions. For example, in the representation where 
H C ont ~ e ijk E^E^F abk is quantized by H cont ~ e ijk F abk (A) d 2 /dA^dAl, solutions are given 
by wave functions ip(A, 7) := TrPexp^ A a (j(t))j a (t) dt, whenever 7 is a smooth, non- 
selfintersecting loop [15]. This happens because the derivatives d/dA bring down two factors 
of the tangent vector, 7 a 7 b , which vanish when multiplied by the antisymmetric tensor F ab . 

It turns out that there exist analogous solutions in the lattice formulation. Take any 
straight "Polyakov loop" a, i.e. a loop without corners that winds around the lattice once, 
and is therefore non-contractible. On a iV 3 -lattice this is a loop a made up of N consecutive 
links in a given direction a, b or c (Fig.l). The contribution of H^{n) to the Hamiltonian 
acting on a wave function ~ Tr V a vanishes at any given site n crossed by a, since it would 
need a wave function with support in at least two independent lattice directions to be non- 
zero. 
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There are 2>N 2 such Polyakov loops (which we will take to be positively oriented), N 2 
in each direction. They will be denoted by a subscript 'p'. The Polyakov loops can be 
parametrized by three integers, namely the coordinates of the point where they intersect one 
of the three planes a = 0, 6 = or c = 0. (Note that each lattice link is contained in exactly 
one Polyakov loop.) We denote the corresponding wave functions by <^(a p ) = TrV^. It is 
easy to see that any wave function ^(cc™) of a multiple Polyakov loop a" = a o a . . . o a 
with winding number n is also annihilated by H . The solution space to the Wheeler- 
DeWitt equation is therefore infinite-dimensional. However, because of the non-linearity of 
the Hamiltonian, it is in general not true that </>(a™)0(/3™) vanishes too. This only occurs 
when (a p ,(3 p ) are a pair of non-intersecting Polyakov loops. We have therefore found: any 
wave function that is a linear combination of terms in 



{(j)(ap 1 )4>((3p 2 ) . . . (f){ojp k ), Op, P p , . . .Up k non — intersecting Polyakov loops} (3.18) 



is annihilated by the discretized Wheeler-DeWitt operator ij . The space of such functions 
is a rather large linear subspace of the original holomorphic Hilbert space. Characteristically, 
its elements are all highly non-local wave functions on the lattice. Note that we need not 
consider separately products of the form (f){a p n )(f){a^) for multiples of the same Polyakov loop 
a p , since these can always be re-expressed via the trace identity (3.10) as sums of elements 
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of (3.18). Whether there exist solutions that are not of this form remains to be explored (see 
also the discussion in the next section). 

From the point of view of non-perturbative quantum gravity, one is interested in the 
structure of the solution space to H^tp® = 0, in particular its scalar product and natural 
self-adjoint operators acting on it. For the subsector of Polyakov wave functions described 
above, there is an induced scalar product from the original Hilbert space of holomorphic wave 
functions. In deriving this inner product, an ambiguity arises because the scalar product on 
the wave functions obtained through the holomorphic transform (2.10) is by construction 
invariant under right and left multiplication by S*?7(2)-matrices, but not bi-invariant under 
SL(2, (D). However, on the explicitly (SL(2, <D)-)gauge-invariant space of Polyakov wave 
functions there is only a small remnant of this gauge covariance: it turns out that the norm 
of a complex wave function depends on the number of link variables that appear in the 
coordinate expression for 4>{a p ). 

This is to be contrasted with the pure SU(2)-case, say. There, in order to simplify 
the computation of scalar products, one often uses a gauge- fixing for a maximal number of 
link variables, which therefore do not any more appear in the calculation. The right- and 
left-invariance of the Haar measure ensures that no physical quantities are affected by this 
choice (see, for example, [16]). In the present case of SL(2,<C), one may also introduce a 
gauge-fixing for some of the links, but one has to keep track of it. Different gauge-fixings 
result in a rescaling of the wave functions. The norm of a Polyakov wave function (f)(a p ) on 
a iV 3 -lattice without any gauge-fixing can be computed using (2.11,12), and is found to be 

\\Ha P )\\sL { 2,v) = e 3Nt/8 . (3.19) 

If m of the N links occurring in a p are gauge-fixed, the norm changes accordingly to 
e 3(N-m)t/8 _ g ca i ar products between (multiple) Polyakov wave functions <fi(a™) and 4>(P P ) 
vanish whenever the underlying Polyakov loops a p and (3 P are different. Furthermore one 
finds that for fixed Polyakov loop a p , different multiples <^(a™) and (j)(a p ) generally have a 
non-zero scalar product, but do not form an orthogonal set, i.e. < <j>{a™), 4>{ot p ) 5 mn . 
Their norms and scalar products can easily be computed using (2.11) and (2.12). 

The next step in the investigation is the search for selfadjoint operators, acting on the 
Hilbert space of the Polyakov wave functions. Natural candidates are the holomorphic trans- 
forms of selfadjoint operators in the S'C/(2)-representation that map Polyakov wave functions 
into themselves. Of course one has to make sure that the final physical expectation values 
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do not depend on a particular gauge-fixing or on the auxiliary parameter t. Our lattice for- 
mulation does for the first time permit rigorous questions about the physical observables and 
their operator spectra. A detailed study of these issues will appear elsewhere. 



4 Gravity on the lxlx 1-lattice 

To illustrate the technicalities involved in the search for solutions other than the Polyakov 
wave functions of the previous section, we now turn to the case of the lxlx 1-lattice with 
periodic boundary conditions. There is a single site s and three (oriented) links emanating 
from it, which we call a, [3 and 7 (Fig.2), to be identified with the 1-, 2- and 3-directions on 
the lattice. 



S . 



7 




We first discuss the gauge-invariant Hilbert space for the 5C/(2)-theory and then apply 
the holomorphic transform (2.9,10). By virtue of the boundary conditions, a, (3 and 7 are 
themselves closed loops. The corresponding holonomy matrices are parametrized according to 
(2.1) by real parameters c^, and ji respectively. The gauge transformations take values in 
a single copy of the gauge group SU (2), located at the site s. For this case we have complete 
control over the gauge-invariant functions, i.e. we can give a complete, non-redundant basis 
of the physical Hilbert space. Following [18], a good set of local coordinates on the six- 
dimensional physical configuration space (SU(2) a x SU(2)p x SU{2) 1 )/ SU(2) S is given by 
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L l («) = ^ Tr ^ = °o 
L 1 (/3) = iTr^ = /? 

L i(7) = J Tr ^r = 7o 

2 i (4-1) 
L 2 (a,(3) = -(Trt^" 1 - TrV a Vp) = a ■ (3 

L 2 (a, 7 ) = j(Tt V a V~ x - Tr V Q V 7 ) = 0-7 

L 2 (/3, 7) = i(TV ^y- 1 - Tr V p V y ) = P ■ 7, 

where we have used a vectorial notation for the "spatial" parameters of the holonomy ma- 
trices, a := (01,02,03) etc. In order to obtain a good global parametrization, one needs to 
add a discrete parameter, which we take to be 



sign(L 3 (a,/3,7)) = sign(a x ^-7). (4.2) 



Also L3 may be written as a linear combination of traces of holonomies involving V a , Vp 
and V 1 [17]. The need for such discrete variables has been repeatedly emphasized by Watson 
[18]. Roughly speaking, the gauge-invariant Hilbert space of square-integrable functions, 
L 2 (SU(2) a x SU(2)/3 x SU(2) 1 )/SU(2) s ,Tr(U l dg)) of the quantum theory is spanned by 
polynomials in the six quantities (4.1) and the discrete variable by a copy of For making 
contact with the corresponding SL(2, (D)-representation, we have to form a Chebyshev basis 
from these polynomials, which is orthogonal and has a simple transformation behaviour under 
(2.10). 

For simplicity, we will first describe the subtheory obtained by reduction to the 1-2- 
plane, say. The relevant classical variables are then Li(a), Li({3) and £2(0,/?). Since the 
real polynomials a^ 1 /3g 2 (a • (3) 113 , regarded as elements of the physical Hilbert space, do not 
have a straightforward holomorphic transform, we use formula (2.11) to identify appropriate 
linear combinations of such polynomials, in order to establish the explicit relation between 
the SU{2)- and the holomorphic SL(2, (D)-representation. For those, one finds 
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[—+11 

* (-l)jo-i ni i( ni _j + ra3 + i)! 
eig [ni ,n 2 ,n 3 ) .- ^ (io _ 1)!(ni + n3)!(ni _ 2jo + 2)! 



jo 

m +1] 



2 J » » 



,_l)fco-i n2 !( n2 _ feo + n3 + i)! 
2. 2 2( fc0 -i) (A;o _i) !(n2 + n 3) ! ( n2 _2A;o + 2)! 

f E E(-!) fel (A) a-- 2 - +2+2fcl 6-- 2fc o™(5^)-- 2i 

\fc 1= 0«i=0 ^ ' ^ ' 

2i-l i-2-1 j i-k-1 ^ 

J\(n 3 -P) g 2(n 2 -fc +n 3 + l-(?) Q 2 (m - j + n 3 + 1 - r) 



(i'i-i) n (^-■-"+2»'* + " +i ) ) 2-(/-"oip-t)i - 

where n, > 0. To give a few simple examples, one has 

eig(0,0,0) = 1 
eig(l,0,0) = a 

eig(2,0,0)=a 2 -^ 

eig(0, 0,1) = a • b 

^(0,0,2) = (a.6) 2 + | + f - ^ 



ezp(l, 1, 1) = a 6oa • 6 
ez£?(2, 1,1)= ao&o a-b-a-b^- 



(4.3) 



(4.4) 



Again, by construction, the highest-order polynomial occurring in eig(ni,ri2,n3) is QJq 1 /?^ 2 (a- 
/3) n3 , and the remaining lower-order polynomials in the sum ensure a simple transformation 
behaviour, which in this case is given by 



C t (eig(n u n 2 ,n 3 )) = e - (ni+n3)(ni+n3+2) * /8 e - (n2+n3)(n2+n3+2) * /8 ei(7(n 1 , n 2 , n 3 ) c . (4.5) 

It is computationally much simpler to determine the action of the Hamiltonian on 
the non-orthogonal basis of holomorphic, square-integrable wave functions {e(ni, ri2, ^3)^ = 
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(chq ) rai ( / 3^) n2 (a c -/3 <D ) n3 }, look for zero-eigenvalue solutions and then use (4.3,5) for the com- 
putation of scalar products, rather than apply on the functions eig(ni,n2,ns)^ directly. 
The quantum Hamiltonian for the 1 x 1 x 1 -lattice consists of a single contribution H®(n) 
and is 

£C = J2 e ijk Tr(V(P &i )T k )Un,*)Pj(n,b)- (4-6) 

a<b 

The plaquette holonomies V(P-r) are defined using the appropriate identifications of links, 
for example, V{P^) = V^V^Vr 1 . 

We now describe our method for finding solutions to the equation Yin m (n) e(n) € = 
for the special case that the wave functions have support only in the 1-2-plane. We first 
calculated the action of the Hamiltonian on a basic holomorphic polynomial e(ni, n 2 , n^)®, 
and obtained 

H € e(ni, n 2 , n 3 ) c = n\n 2 ( e{n\ - 1, n 2 - 1, ra 3 + 3) c - e{n\ — 1, n 2 — 1, n 3 + 1) € ) 
+(nin 2 + nin 3 + n 2 n 3 - n 3 ) ( e(n 1 ,n 2 , n 3 ) c - e(ni, n 2 , n 3 + 2) c ) 
+(nin 2 - nin 3 ) e(ni - 1, n 2 + 1, n 3 + l) c + (nin 2 - n 2 n 3 -) e(ni + 1, n 2 — 1, n 3 + l) c 
+(nin 2 + nin 3 + n 2 n 3 - n 3 ) ( e(ni + 2, n 2 + 2, n 3 ) c - e(ni + 2, n 2 , n 3 ) c 
-e(n 1; n 2 + 2, n 3 ) c ) - (nin 2 - nin 3 - n 2 n 3 — n§ — 2n 3 ) e(ni + 1, n 2 + 1, n 3 + 1) € 
+nin 3 ( e(ni - 1, n 2 + 1, n 3 - l) c - e(ni - 1, n 2 + 3, n 3 - 1) € ) 
+n 2 n 3 ( e(ni + 1, n 2 - 1, n 3 - I) 10 - e(n x + 3, n 2 - 1, n 3 - 1) € ) 

-(2nin 3 + 2n 2 n 3 + n\ + 2n 3 ) e(ni + 1, n 2 + 1, n 3 - l) c + (nin 3 + 2n 2 n 3 + nj + 2ra 3 ) 

e(ni + 3, n 2 + 1, ra 3 - 1) € + (2nin 3 + n 2 n 3 + n§ + 2n 3 ) e(rti + 1, n 2 + 3, n 3 - l) c 
-(nin 3 + n 2 n 3 + n| + 2n 3 ) e(ni + 3, n 2 + 3, n 3 - l) c 

-n 3 (n 3 - 1) ( e(m + 2, n 2 , n 3 - 2) c - e(ni + 4, n 2 , n 3 - 2) c + e(m, n 2 + 2, n 3 - 2) c 
-3 e(m + 2, n 2 + 2, n 3 - 2) c + 2 e(ni + 4, n 2 + 2, n 3 - 2) c - e(m, n 2 + 4, n 3 - 2) c 
+2 e(ni + 2, n 2 + 4, n 3 - 2) c - e(ni + 4, n 2 + 4, n 3 - 2) c ), 

(4.7) 

for rii > 0. That is, independent of the values of the nj, the right-hand side of (4.7) is a linear 
combination of (at most) 26 terms, and the arguments of the e{n\, n 2 , n 3 ) c occurring lie in 
the ranges —1 < Ani < 4, — 1 < An 2 < 4 and —2 < An 3 < 3. 
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Since the functions e(ni,n 2 ,n 3 ) € form a (non-orthogonal) basis for the gauge-invariant 
Hilbert space, we can reformulate the zero-eigenvalue condition as an infinite set of equations 
for the m(rei, n 2 , n 3 ) obtained by setting the coefficient of each e(n\, n.2, n 3 ) on the right-hand 
side of Yin- H € m(ni,n2,n^) e(ni,n 2 ,n 3 ) c to zero. The general condition can be labelled by 
three integers and is easily derived from (4.7), yielding 

C[ni, n 2 , n 3 ] := (m + l)(n 2 + 1) ( m(n\ + 1, n 2 + 1, n 3 — 3) - m{n\ + 1, n 2 + 1, n 3 — 1) ) 

- (n 2 n 3 - nin 2 - n\ + n 3 ) m(ni - 1, n 2 + 1, n 3 - 1) + (nin 2 + nin 3 + n 2 n 3 - n 3 ) 
m(ni, n 2 , n 3 ) - ((m - 2)n 2 + (ni - 2)n 3 + n 2 n 3 - n§) m(ni - 2, n 2 , n 3 ) 

- (nin 2 + ni(n 3 - 2) + n 2 (n 3 - 2) — n 3 + 2) m(ni, n 2 , n 3 - 2) 

- (— nin 2 + nin 3 - n 2 + n 3 ) m(ni + 1, n 2 — 1, n 3 - 1) 

- (m (n 2 - 2) + nin 3 + (n 2 - 2)n 3 — n 3 ) m(ni, n 2 — 2, n 3 ) 

- {n\U2 - n\n 3 - n 2 n 3 - n\ + 2n 3 ) m(ni - 1, n 2 — 1, n 3 - 1) 

+ (ni + l)(n 3 + 1) ( m[n\ + 1, n 2 - 1, n 3 + 1) - m(ni + 1, n 2 — 3, n 3 + 1) ) 
+ ((m - 2)(n 2 - 2) + (m - 2)n 3 + (n 2 - 2)n 3 - n\) m(rai - 2, n 2 - 2, n 3 ) 

- (2nin 3 + 2n 2 n 3 + nl + 2n\ + 2n 2 - 1) m(n\ - 1, n 2 — 1, n 3 + 1) 

+ (nin 3 + 2n 2 n 3 + n\ + ri\ + 2n 2 - n 3 — 2) m(ni - 3, n 2 — 1, n 3 + 1) 
+ (2riin 3 + n 2 n 3 + n\ + 2ni + n 2 - n 3 - 2) m(ni - 1, n 2 - 3, n 3 + 1) 

- (nin 3 + n 2 n 3 + n| + m + n 2 - 2n 3 - 3) m[n\ - 3, n 2 — 3, n 3 + 1) 

+ (n 2 + l)(n 3 + 1) (m(ni - l,n 2 + l,n 3 + 1) - m(n\ - 3,n 2 + l,n 3 + 1) ) 

- (n 3 + 3n 3 + 2) ( m{n\ - 2, n 2 , n 3 + 2) - m[n\ - 4, n 2 , n 3 + 2) + m(ni, n 2 - 2, ra 3 + 2) 

- 3 m(ni - 2, n 2 - 2, n 3 + 2) + 2 m(ni - 4, n 2 - 2, n 3 + 2) - m(n\ , n 2 - 4, n 3 + 2) 
+ 2 m(ni - 2, n 2 - 4, n 3 + 2) - m(ni - 4, n 2 - 4, n 3 + 2) ) = 0. 

(4.8) 

Notice first that the two-dimensional wave functions e(ni,n 2 ,n 3 ) fall into two sectors 
which are mapped into themselves under the action of the Hamiltonian. These are (i) the 
even sector: either all n, are even or all are odd; (ii) the odd sector: two rii are even and 
one Hi is odd or vice versa. Hence it suffices to investigate the two sectors separately. 

The special solutions discussed in the previous section correspond to all functions of 
the form e(ni,0,0) (E and e(0, n 2 ,0) (C , and (4.8) imposes no conditions on the corresponding 
coefficients m(ni,0, 0) and m(0,re 2 ,0). We will in the following set these coefficients to zero, 
because we are interested in the possible existence of other solutions. In order to tackle the 
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system of equations in a well-defined manner, it is useful to define the order ord(C[ni, ri2, n^]) 
by 



ord(C[n 1 ,n 2 ,n 3 \) := max{ni + n 3 , n 2 + n 3 }, (4.9) 

and then try to solve simultaneously, order by order, the sets of equations (4.8) of the same 
order, eliminating coefficients m(ni, 712,^3) of higher order (with 'order' defined analogously 
to (4.9)) in terms of the lower order ones. We have investigated the even sector and solved 
iteratively in this manner the order-0, 2, 4, 6, 8, 10 and 12 equations. There are 1 order-0 
equation, 5 order-2, 13 order-4, 25 order-6, 41 order-8, 61 order-10 and 85 order-12 equations. 
Taking into account equations up to this order, one finds that for any solution to the Wheeler- 
DeWitt equation the coefficients have to satisfy simultaneously 



order 2 : m(0, 0, 2) = m(l, 1, 1) = m(2, 2, 0) = 

order 4 : m(l, 3, 1) = m(3, 1, 1) = m(l, 1, 3) = m(2, 0, 2) = m(0, 2, 2) = m(0, 0, 4) = 

1 38 

m(2, 4, 0) = m(4, 2, 0) = - m(2, 2, 2), m(3, 3, 1) = — m(2, 2, 2), 

no conditions on m(2, 2, 2), m(4, 4, 0) 

7 2 
order 6 : m(l, 5, 1) = m(5, 1, 1) = -- m(2, 2, 2), m(l, 1, 5) = — m(2, 2, 2), 

m(4, 0, 2) = m(0, 4, 2) = m(2, 2, 2), m(3, 1, 3) = m(l, 3, 3) = ^ m(2, 2, 2), 

m(0,0,6) = -^L m (2,2,2), 

(4.10) 

i.e. the coefficients m of order and 2 are completely fixed, but at order 4 there appear two 
free parameters, m(2, 2,2) and m(4, 4, 0). (Of course we could have solved in terms of other 
order-4 parameters.) Unfortunately we cannot be sure whether the conditions C[n±, ri2, 713] = 
for ord(C) > 12 do not (through coupling among equations of different order) fix these 
parameters, although the behaviour of the conditions evaluated so far does not make it seem 
likely. We conjecture that there is an infinite number of free parameters (with an increasing 
number of free parameters at each order), corresponding to an infinite set of solutions to the 
discretized Wheeler-DeWitt equation (beyond those coming from the Polyakov loops) . These 
would be of the form of (presumably infinite) linear combinations, parametrized by those free 
parameters. Clearly then the question arises of whether these solutions have a finite norm. 
Since we have not even found a single explicit solution of this type, we are unable to answer 
this question presently. 
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We will now have a brief look at the full 1 x 1 x 1-lattice theory. However, we will 
not attempt to solve the eigenvalue equation directly, since the analogues of (4.7) and (4.8) 
contain about ten times as many terms. The discussion is meant to serve as an illustration 
of how to set up the gauge-invariant Hilbert space in a symmetric way, and formulate the 
eigenvalue problem in principle. The problem is the incorporation of the classical discrete 
degree of freedom sign(L3(a, (3, 7)) in the quantum theory. Although the functions 



{ a ^p^(d-pr\d-jr0^r( S i g n(dxp^)y, ^ = 0,1,2,..., * = o,i} (4.11) 

could in principle serve as a basis for the gauge-invariant Hilbert space, it is difficult to 
incorporate sign(a x (3 ■ 7) in integrations, since it is not a smooth function on the classical 
group manifold. Another admissible choice is 

{a^/3 " 2 7o 3 (5 • (3) n *(a ■ ^(a x (3 • tT 6 , ^ = 0, 1, 2, . . .}, (4.12) 

but this is not symmetric with respect to the three lattice directions. A choice that solves 
both of these problems is 

{e(m, n 2 , n 3 , n 4 , n 5 , n 6 , z) := a^/^o 3 (a ■ ^ (5 ■ tT 5 ■ iT 6 (5x^ 7 T, 

Hi = 0,1,2,..., z = 0,1}. 

It contains all gauge-invariant information about the original Hilbert space, without 
being overcomplete. There is a corresponding, non-orthogonal basis for the square-integrable, 
holomorphic wave functions, obtained by substituting as usual the real parameters in (4.13) 
by their complex counterparts. Acting with the Hamiltonian iJ € on such a state may produce 
terms that contain higher-order powers (a® x (3® ■ J®) n , n > 1, which then have to be re- 
expressed as a sum of (complex) terms of the form (4.13) using the identity 

(5 C X (3® . f € } 2 = (1 _ (5 € )2)(1 _ ( ^C )2)(1 _ (f € )2) _ (1 _ (5 C)2 )( £<D . f € )2 

- (1 - 03 C ) 2 )(5 C • 7 C ) 2 - (1 - (7 C ) 2 )(5 C • ^? + 2(a € • ^ c )(a c • 7 C )(^ C ■ 7 C )- 

(4.14) 
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This is an unambiguous prescription and leads to a Hamiltonian action that maps states 
of the form m(ni, 112, 113, 114,115, iiq, z)e{yi\,ni, n.3, ns, n§, z)® into themselves. However, 
as already mentioned, the Hamiltonian is a lengthy expression and solving the zero-eigenvalue 
problem is certainly not a straightforward task. This is probably just a reflection of the non- 
triviality of quantum gravity. Still we cannot rule out the existence of a basis for the gauge- 
invariant Hilbert space that leads to a simplification of the eigenvalue problem. However, 
this would presumably be an overcomplete basis, which then leads to problems of a different 
kind. Firstly, one cannot just set the coefficients of each "basis" wave function to zero, and 
secondly one has to eliminate spurious solutions by hand. By contrast, an advantage of our 
formulation is that the general condition on the wave function coefficients m(n) is known 
explicitly, so that one has full control over all the relevant physical parameters. This is 
important if one for instance decides to introduce a cut-off in the Hilbert space, in order to 
approximate the problem of finding zero-eigenvectors by a finite-dimensional one. 



5 Conclusions 

We have described above a regularized version of non-perturbative canonical quantum 
gravity on a cubic lattice with periodic boundary conditions. The quantum Hamiltonian on 
the lattice acts combinatorially on holomorphic wave functions labelled by lattice loops. The 
lattice represents an entire diffeomorphism equivalence class and the formalism is manifestly 
gauge- and diffeomorphism-invariant at the kinematical level. With a specific choice for 
the discretized Hamiltonian and a factor ordering in the quantum theory, we are able to 
identify an infinite-dimensional space of solutions to the discretized Wheeler-DeWitt equation, 
which moreover have finite norm with respect to a natural scalar product on the space of 
holomorphic SL(2, (D)-functions. The solution space is labelled by global Polyakov loops and 
their multiples, and corresponds to ~ 3N 2 physical degrees of freedom (compared to ~ 6iV 3 
before imposing the Hamiltonian constraint). Still, more research is needed to determine 
whether this exhausts the solution space. For the example of the 1 x 1 x 1-lattice, we 
illustrated how one may go about a systematic search for more solutions. Although we have 
not yet been able to find any, preliminary results suggest there may be an infinite set of 
solutions beyond the Polyakov ones. However, even if that is the case, it may still happen 
that they are not square-integrable. 

Since our lattice regularization is rather different from the point-split regularization used 
in the formal continuum approaches, it is remarkable that our solution space is reminis- 
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cent of the smooth, non-intersecting solutions of [6] (thus suggesting that these in fact are 
not "spurious"). We expect this feature to be fairly robust under a change of the classical 
discretized Hamiltonian, because the existence of the solutions depends only on the antisym- 
metric structure of e l i k pjPk- However, we have not found straightforward analogues of the 
intersecting-loop solutions of [15,19]; this can be traced back to the non-locality of the lattice 
Hamiltonian. The great advantage of our regularization is the existence of a well-defined 
scalar product at every stage. This will be crucial in all further investigations of the solution 
space. For example, it would be interesting to understand how a different factor ordering of 
H € chan ges our results. 

In order to avoid confusion, it should be pointed out that we are proceeding somewhat 
differently from the path outlined (for the continuum theory) by Ashtekar and collabora- 
tors [20,13]. They propose to solve the gauge and diffeomorphism constraints within a real 
framework, based on the spin connection r* , and then go to a holomorphic representation 
to solve the Hamiltonian constraint, whereas our formulation takes place entirely within the 
complex formulation based on A\ = T l a — i K l a (where K l a is related to the extrinsic curvature 
via K l a = K a \yE bl ). This does not exclude that a close relation may emerge at the level of 
dynamics. 

Finally, as already mentioned in the introduction, one has to face the question of the 
continuum limit (not to be confused with a weak-field limit) of the regularized theory. The 
diffeomorphism invariance of general relativity makes this a fundamentally different issue 
from that in lattice gauge theory. If one wants to avoid bringing back in an ultraviolet 
cutoff a, the only free parameter is the lattice size N, and one would expect that in the 
limit for growing N a prospective continuum theory is approximated ever better. Since (at 
least part of) the solution space is known for every finite N (and the iV-dependence of our 
construction is rather explicit), one might investigate the limit as N — > oo of these spaces 
directly. This will probably become more meaningful once observables and possibly matter 
fields have been included, so that one can study their spectral properties as a function of 
N. Our construction suggests that in such a "continuum" limit some fundamental discrete 
structure is retained, although (as pointed out in [21]) this does not necessarily preclude the 
appearance of divergences. 
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